T 


Blowing in the Wind: II. Creation and Redistribution of Refractory Inclusions 

in a Turbulent Protoplanetary Nebula 

Jeffrey N. Cuzzi, Sanford S. Davis, and Anthony R. Dobrovolskis 
Space Science Division, Ames Research Center; cuzzi@cosmic.arc.nasa.gov 
Submitted to Icarus, April 14 2003 

Abstract: Ca-Al rich refractory mineral inclusions (CAIs) found at 1-6% mass 
fraction in primitive chondrites appear to be 1-3 million years older than the 
dominant (chondrule) components which were accreted into the same parent 
bodies. A prevalent concern is that it is difficult to retain CAIs for this long 
against gas-drag-induced radial drift into the sun. We reassess the situation in 
terms of a hot inner (turbulent) nebula context for CAI formation, using ana- 
lytical models of nebula evolution and particle diffusion. We show that outward 
radial diffusion in a weakly turbulent nebula can overwhelm inward drift, and 
prevent significant numbers of CAI-size particles from being lost into the sun 
for times on. the ojcder of 10 6 years, C AIs can form early, when the inner, nebula 
was hot, and persist in sufficient abundance to be incorporated into primitive 
planetesimals at a much later time. Small 0.1 mm diameter) CAIs persist 
for longer times than large (i> 5mm diameter ones. To obtain a quantitative 
match to the observed volume fractions of CAIs in chondrites, another pro- 
cess must be allowed for: a substantial enhancement of the inner hot nebula in 
silicate- forming material, which we suggest was caused by rapid inward drift of 
meter-sized objects. This early in nebula history, the drifting rubble would have 
a carbon content probably an order of magnitude larger than even the most 
primitive (Cl) carbonaceous chondrites. Abundant carbon in the evaporating 
material would help keep the nebula oxygen fugacity low, plausibly solar; as in- 
ferred for the formation environment of CAIs. The associated production of a 
larger than canonical amount of C0 2 might also play a role in mass-independent 
fractionation of oxygen isotopes, leaving the gas rich in 16 O as inferred from 
CAIs and other high temperature condensates. 

1 Introduction 

Chondrite parent bodies are dominated by particles with a surprisingly well-defined range 
of physical, chemical, and petrographic properties. Fe-Mg-Si-0 mineral chondrules, many of 
which solidified from melted precursors stable at T < 680K, constitute 30-80% of primitive 
meteorites. Most workers in the field believe that chondrules are formed by either localized or 
nebula scale energetic events operating on freely floating precursors of comparable mass, at 
some location or locations in the protoplanetary nebula (see eg. Grossman 1989, Grossman 
et al. 1989, Boss 1996, Connolly and Love 1999, Jones et al. 2000, and Desch and Connolly 
2002). 

A different class of mineral grains called Ca-Al-rich inclusions (CAIs), whose constituent 
minerals condense out of nebula gas at a much higher temperature ( T > 1500K), make up 
1-6% of primitive meteorites depending on CAI and chondrite type. Many, but not all, have 
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apparently been subsequently melted (Grossman et al. 2002). The formation locale of CAIs 
is also unknown; some place it in the hot inner part of the terrestrial planet zone where the 
meteorites themselves subsequently evolve {eg., Cassen 2001) and others in the immediate 
solar neighborhood {eg., Shu et al. 1996). There is increasingly strong radioisotope evidence 
that most CAIs are several 10 6 years older than most chondrules (Mostefaoui et al. 2001, 
Huss et al. 2001, Amelin et al. 2002). How these high-temperature mineral fragments find 
themselves intimately mixed with lower-temperature minerals, especially after such a long 
time, has been a puzzle. In fact, many have questioned the apparent age differences merely 
on the belief that CAIs would be lost by gas drag into the sun in a time much shorter than 
the apparent age difference. 

In this paper we address the survival of CAIs for times > 10 6 years. We show that radial 
diffusion counteracts the expected inward drift, given a certain amount of (weak) nebula 
turbulence. CAIs might form in the terrestrial planet zone - but early, at a high temperature 
- and resist inward drift processes for more than 10 6 years. CAIs might even form much closer 
to the sun and still diffuse “upstream”, radially outwards into cooler chondrule formation 
and accumulation regions. A combination of the two different source regions is both possible 
and consistent with the properties of the two main types of CAIs. 

As part of this work, we describe a new process by which the inner solar nebula might be- 
come significantly enhanced in rock-forming elements, CO, and CO 2 relative to over nominal 
abundances. The process is that of inward drift of very primitive, carbon-rich meter-sized 
particles into regions which are hot enough to evaporate them. Halley dust and cometary 
IDPs, for instance, have C/O ratios close to solar - that is, enriched by an order of magni- 
tude relative to Cl chondrites (Jessberger et al. 1988, Jessberger and Kissel 1991; Lawler 
and Brownlee 1992). Evaporation of this material along with ferromagnesian silicates allows 
the silicate abundance to increase without elevating the oxygen fugacity, as would occur for 
a nominal Cl mix. Furthermore, Thiemens (1996, 1999) has argued that isotope exchange 
reactions involving C0 2 fractionate heavy oxygen isotopes into the C0 2 , leaving reactive 
oxygen enriched in 16 O. The C0 2 is non-reactive at these temperatures, and carries the 
heavy isotopes away, leaving high-temperature condensates to become enriched in 16 0 as 
observed. Quantitative development of these chemical and isotopic roles of carbon are be- 
yond the expertise of these authors, so these aspects are not developed further but merely 
mentioned as plausibility arguments for others to pursue. In the remainder of this section we 
mention some basic aspects of turbulence, and of the influence of gas drag on particles. In 
section 2 we describe our model which includes a globally evolving nebula, an inner nebula 
“CAI factory”, and subsequent radial diffusion. In section 3 we present our model results. 
In section 4 we summarize and discuss the implications. 

1.1 Turbulence 

Weak global turbulence is an essential aspect of this scenario. It remains a subject of debate 
whether the nebula gas was turbulent or laminar during the CAI-chondrule era (Stone et 
al. 2000, Richard and Zahn 2001, Klahr and Bodenheimer 2002, Richard 2003). We have 
suggested that the typical chondrule size and size distribution themselves may be indicators 
of weak nebula turbulence (Cuzzi et al. 2001). In that paper, we also pointed out that 
adequate turbulence for interesting effects on particle evolution might be a byproduct of 
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nebula evolution even if it is not a driver, and that only about 10~ 5 of the gravitational energy 
which is continually released as the nebula evolves needs to be tapped to generate this amount 
of turbulence. Turbulence can be described by its kinetic energy per unit mass V g /2, and the 
sizes of its largest, or integral scale eddy L and its smallest, or Kolmogorov scale eddy rj. The 
turbulent Reynolds number Re = ( L/r ?) 3 / 4 , is also written as Re = vt/v = V g L/u = acH/u, 
where ut and v are the turbulent and molecular viscosities and H is the gas vertical scale 
height. V g is the velocity of the large eddies, which contain most of the energy; the turnover 
time of these eddies, ti = 1/SIl — L/V g , is generally assumed to be the orbit period. We 
will not distinguish here between the gas turbulent diffusivity V and turbulent viscosity 
ut (Prinn 1990, Stevenson 1990; cf. also Cuzzi et al. 2001 for more discussion). The 
smallest eddies have the smallest velocities and the shortest turnover times (Cuzzi et al. 
2001, Cuzzi and Hogan 2003). In reality, turbulence was probably radially, vertically, and 
temporally variable, but those specifics are beyond the scope of this paper, if not of current 
understanding; here we adopt two very simple, analytical models of turbulent intensity and 
viscosity which we expect will span a sufficiently wide range of behavior to provide insight 
into, what , may have happened in the actual nebula (discussed in section 2). First,. ,we sketch 
below the particle-gas dynamics involved. 

1.2 Particle drift and diffusion 

Particles such as CAIs and chondrules, which are smaller than the gas molecular mean free 
path, have a gas drag stopping time 

t s = r pPs /cp g , (1) 

where r p is particle radius, p s is particle material density, c is the nebula sound speed, 
and p g is the nebula gas density (Weidenschilling 1977). Particles equilibrate with gas 
velocities on the timescale t s . For CAIs and chondrules, and canonical nebula parameters, 
t s /ti = St « 1. The value of St (the Stokes number) entirely determines the interaction 
of the particles with a turbulent gas (Volk et al. 1980, Markiewicz et al. 1991, Cuzzi et al. 
1993, Cuzzi and Hogan 2003). 

Weidenschilling (1977) presented radial drift velocities and lifetimes for various size par- 
ticles in a number of nonturbulent nebula models, neglecting collective effects of adjacent 
particles (see however Nakagawa et al. 1986 and Cuzzi et al. 1993). The angular velocities 
of small particles (St << 1) are forced to nearly match that of the gas, which orbits more 
slowly than Keplerian because it experiences a net outward pressure gradient acceleration 
A g — - p dP/dR. Thus, the particles experience an uncompensated inward acceleration A g 
and acquire an inward drift velocity of 

1 Up T/2 

v d = A gt, — t, « -27— f., (2) 

p dr r 

where 7 is the ratio of the pressure gradient to the central force gravity and V K is the 
Keplerian velocity. It turns out that 7 w 2 x 10“ 3 is valid for a wide range of plausible 
nebula locations and conditions (Cuzzi et al. 1993). Then 

v d « ~ - 2 jV K t s Cl « -27 V K St ( 3 ) 

r r\l 
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in the regime of interest where pd « p g (cf. Nakagawa et al. 1986, equation 2.11). Thus, 
for small particles, V d increases with St, until St ~ 1 . The most rapidly drifting particles are 
those with St = 1; these drift at V d — ~-^Vk {eg-, Weidenschilling 1977). It happens that unit 
density, meter-radius particles satisfy St ~ 1 nearly throughout the entire nebula, because 
variations in Vt and drag coefficient offset variations in p g : We refer to the drift velocity of 
“large” meter-sized particles as V L — —yV K . For very large particles with t 3 Cl = St » 1, 
the drift velocity again decreases, as V d = —27 V^/St. 

The second effect of relevance is that the gas within which the particles are drifting is 
itself evolving radially due to angular momentum transfer. Lynden-Bell and Pringle (1974; 
henceforth LBP) showed that inner regions evolve inwards (into the sun) and outer regions 
evolve outwards (conserving angular momentum), with the boundary between these regions 
itself evolving outwards with time. In our actual models (section 2), we will use the analytical, 
radially dependent results of LBP, but for background purposes in this section the radial 
drift of the nebula gas V n may be approximated by an inward velocity 


V n = 


o t'7’ 
2 ~ 


3 udl 

_____ 
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(Lin and Papaloizou 1985, Cassen 1996, 2001; Gail 2001). For nebula parameters, both V d 
and V n are in the 10 cm sec -1 range (discussed in more detail below). 

The third effect (if turbulence is present), is radial diffusion of the particles. The random 
inertial space velocity of the particles V p determines the particle diffusion coefficient V = 
1 'rVp/Vg', Volk et al. (1980) and Cuzzi and Hogan (2003) show that V p = V p /(1 + St) « 
to a very good approximation for CAIs and chondrules, so V = v T /( 1 4- St) « i/ T (see also 
Cuzzi et al. 1993, Appendix B). That is, because the stopping times of mm-sized particles 
are so much shorter than the turnover times of the most dispersive eddies, the velocity and 
mixing lengthscales of these particles are comparable to those of gas molecules. 

A quick comparison of these three characteristic velocities reveals that V p is much larger 
than either V n or V d for mm-to-cm-sized particles in the regime of interest. Specifically, 
for these particles V p « V g = ca V 2 {eg., Cuzzi et al. 2001). For nebula a in the range 
of 1CT 4 to 10 -2 , the random, turbulence-driven particle velocities V p are much larger than 
drift velocities due to either gas drag or inward nebula evolution ( V p /V d and V p /V n both 
~ 10 2 - 10 4 ). 

The main point is the radial mass flux of particles under either diffusion or advection. 
Suppose the vertically averaged CAI concentration C decreases outwards, with a radial 
lengthscale A, as expected if CAIs form in some hot inner region. Then the ratio of (outward) 
diffusive flux to (inward) advective mass flux is 


Va{dC / dr) Va{C/\) r 

{V n + V d )aC ~ V n aC ~ A 


where we have made use of the fact that V d « V n for mm-sized particles, and have substituted 
V n reVjr (equation 4). 

Equation 5 shows that diffusive (outward) mass flux dominates advective (inward) drift 
until the radial scale of variation of the concentration is larger than the radius itself - essen- 
tially, a nearly constant radial concentration profile. It is no more meaningful to discuss radial 
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evolution of such particles without considering radial diffusion than it is to discuss “settling 
to the midplane” of such .small particles without considering vertical diffusion (Dubrulle et 
al. 1995, Cuzzi et al. 1993, 1996). Indeed, Weidenschilling (1977) mentions some caveats on 
this subject, and it has also been discussed in detail by Morfill and Volk (1984), but with 
an approach that suppresses the radial transport aspects we wish to emphasize. In related 
research, Bockelee-Moravan et al. (2002) used a full numerical diffusion-advection model to 
show that tiny grains from the inner solar system, which do not incur appreciable gas-drag 
radial drift, diffuse throughout the solar system on short timescales. 


2 Method of solution 


As a first step towards a solution of this complex problem, we frame a scenario for the “CAI 
factory” and adapt analytic solutions already in the literature to assess how the distribu- 
tion of these particles evolves radially with time. We first describe our analytical nebula 
models (section 2.1) and our scenario for the CAI factory, including the new aspect of drift 
enhancement of the inner regions in silicates (section 2.2). We then describe how we treat 
simultaneous drift and diffusion (section 2.3). We present typical results of the model in 
section 3. 


2.1 Nebula Evolution 


We must allow for the fact that, over the few x 10 6 year duration of interesting evolutions, the 
mass density and radial velocity of the nebula will evolve significantly by virtue of angular 
momentum and mass redistribution. We use the LBP analytic solutions for the radial profile 
of the local surface mass density cr g (r,t ) and the nebula radial evolution velocity V n (r,t) in 
an evolving nebula. With the exception of the short-lived, hot inner region, we assume that 
silicates represent a constant fraction of the gas surface density. 

The LBP formalism allows for a variety of powerlaw radial profiles of ux = acH = 
ac 2 /0 (see also Hartmann et al. 1998). We choose two limiting cases; case 1 is a constant 
viscosity case vxi = ac 2 /Q.d and case 2 only has constant ac 2 , which results in u T 2 = 
ac 2 /S7c(r/r£)) 3 / 2 oc r 3 / 2 . fl# is the orbit frequency at a reference disk radius rp defined 
below. These two cases bracket a range of plausible r— dependencies of v ■?. R. Bell (personal 
communication, 2002) has found that, in her nebula energetics models (Bell et al. 1997), 
the radial dependence of visosity is actually not too far from one or the other of these 
dependencies (over different radial regions, the boundary between which changes with time). 
The LBP solutions may be exactly rewritten in simpler forms. For constant up: 

(M d \( t \~ 5 ' 4 [ ( r/r D ) 2 

" i(r ' i)= te)( l+ d exp r (1^ 


( 6 ) 


and for constant ac 2 : 
a 2 (r, t ) 


KiM) = “r 


r / 2 K) 


i ~ 


{r/r D )‘ 


(1 4- t/t e ) 



-2 


exp 


(r/r g ) 1/2 
(1 -f t/t e ) 


(7) 
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Figure 1: Surface mass density of the nebula gas, for the two different models we will study 
in this paper: constant viscosity models (solid lines) and constant ac 2 models (dashed lines). 
The mass density decreases with time; the curves begin at our first timestep (7 x 10 4 years), 
and then are plotted on even multiples of 0.5 Myr until- 3 Myr years (3 Myr curves are plotted 
with a heavier weight). This set of evolutions is for nebula a = 3x 10 4 . The case 2 viscosity 
is multiplied by a factor of 3 to keep the mass fluxes of the two models (figure 2) in good 
agreement at each timestep. 


V n2 (r,t) 


3 Mrp) 

2r 



(r/r g ) 1/2 ~ 
(1 4- t/t e ). 
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In the above equations, Mp is the initial disk mass, M 0 the solar mass, and G is the 
gravitational constant. The reference radius rp is obtained using the angular momentum per 
unit mass scaling of Cassen (1996): rp = hl/GM Q , with h 0 the specific angular momentum 
of the solar system. The disk evolutionary timescale is t e = r^Qp/bac 2 , with 6=12 for 
the constant ut case and b = 3/4 for the constant ac 2 case. These equations may be of 
general use, since we have converted the arbitrary mathematical constants of LBP into three 
meaningful physical properties of the disk: Mp, rp, and a. Specifically, we make the following 
substitutions for LBP’s a and C: a = (GM Q rp)~ 2 (relating radii to the nebula’s specific 
angular momentum); and Cjv T = 3Mp/r 2 D (from integrating over the disk mass). A more 
complicated analytical model, although essentially based on the same three parameters, has 
been developed by Stepinski (1998). Davis (2003) also describes several different approaches 
to analytic modeling. 

Typical results from such models are shown in figures 1-4, and model parameters are 
given in Table 1. In all figures, the constant viscosity models are shown in solid lines and 
the constant ac 2 models in dashed lines. In spite of the significant differences in the radial 
a g profiles, the radial profiles of mass influx are quite similar and are nearly independent of 
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Figure 2: Radial mass flux of material (absolute value) as a function of radius and time, for 
the two models of figure 1 . The mass flux decreases with time; values are plotted at the 
initial timestep (7 x 10 4 years) and then every 0.5 Myr from 0.5 Myr until 3 Myr (heavy 
curves at 3 Myr). Inwards of the cusp the mass flux is negative (inwards). Inwards of 
the cusp, mass flux is inwards (negative); outwards of the cusp, the mass flux is positive 
(outwards). The mass fluxes for the two cases are aligned by a factor of 3 increase in the 
viscosity of the constant ac 2 model relative to a = 3 x 10 -4 (section 2.1.) 
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Figure 3: Nebula gas velocities as functions of radius and time, for the two viscosity models 
(solid: constant viscosity models; dashed: constant ac 2 models). Upper plot: gas drift 
velocities for a = 3 x 10 -4 , lower plot: gas drift velocities for a = 3 x 10 -3 






r (AU) 

Figure 4: Millimeter-radius (unit density) particle drift velocities as functions of radius and 
time, for the two viscosity models (solid: constant viscosity models; dashed: constant ac 2 
models). Upper plot: particle drift velocities for a = 3 x 10 -4 , lower plot: particle drift 
velocities for a = 3 x 10 -3 . The larger drift velocities are caused by the fact that the nebula 
evolves more rapidly at higher a values, and thus the gas density is lower at any given time, 
resulting in a larger t s and a larger (see equation 2). 
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radius in the inner solar system. The similarity in mass flux profiles (figure 2) is related 
to the different radial velocity profiles (figure 3). In fact, it is a natural feature of the 
LBP solutions that the radial dependencies of a g and v? nearly cancel throughout the mass 
inflow region. The nearly constant mass influx depends on their product (cf., eg., Lin and 
Papaloizou 1985). This suggests, for instance, that <? g (r,t ) would be proportional to r _1 for 
a viscosity vt oc r such as suggested by Hartmann et al. (1998) or proportional to r -1 / 2 
for the shear-driven /3 - viscosity profile vt oc r 1 / 2 suggested by Richard and Zahn (1999) (see 
also Davis 2003). The mass flux and overall nebula mass decrease with time as expected, 
with magnitudes (in solar masses per year) which are typical for nebulae of age 10 5 — 10 6 
years (radially outwards of the cusps in figure 2, the mass flux reverses and is outwards.) 

2.2 The CAI factory 

For the purpose of this modeling exercise, we adopt a highly simplified picture of the en- 
vironment in which CAIs may have formed and evolved. Such a simplified model can only 
be asked to address the crudest possible characterization of the OAls themselves. (JAls are 
primarily oxides and silicates of calcium, aluminum, and magnesium, with small amounts 
of titanium and other refractory elements (MacPherson et al. 1989; Brearley and Jones 
2000). Some CAIs appear to be, essentially, condensates or aggregates of condensates from a 
cooling high temperature gas, chemically modified to greater or lesser degree by subsequent 
reaction with the surrounding gas as it continues to cool, and/or subsequent but limited 
duration heating events (eg., MacPherson and Davis 1994, especially their figure 11). Other 
kinds of CAIs are “igneous” in that they seem to have solidified from a melt - possibly a 
melt of earlier condensates (eg., MacPherson et al. 1988, Grossman et al. 2002). The most 
widely occurring type of CAIs (type A), which are found in essentially all meteorite classes, 
are the smallest, contain the highest temperature minerals', and have the closest structural 
appearance to simple condensates; the most abundant (type B) CAIs, which are found al- 
most entirely in one (CV) chondrite class, are larger, tend to be more dominated by lower 
temperature (but still relatively refractory, low-Silicon) minerals, and tend more to manifest 
recrystallization from a melt than condensation from a gas. 1 

The condensation- reaction process has not continued to completion in either case (G. 
MacPherson, personal communication, 2002), arguing against a simplistic formation process 
in which CAIs are always in equilibrium with the same parcel of gas as it cools over infinite 
time. It will ultimately be necessary to relate timescales of nebula radial evolution (outwards, 
into cooler regions) to timescales of chemical equilibration of the object with the surrounding 
gas. Our model offers the potential for quantitative assessments of this type. 

While complications are easy to imagine, (multiple, relatively brief heating events, vertical 
thermal structure, and so on), we present a simple picture based on prior models of inner 
nebula energetics and temperature as an initial context for pondering these two very different 
refractory particle types. Our simple picture (figure 5) is that the innermost nebula is 
sufficiently hot at early times out to some boundary for all silicates and refractory oxides 

Vigorously, “type A” and “type B” refer to specific mineralogies, traditionally those found in CV chon- 
drites, as well as to their characteristic size, porosity, and petrography. Here we will suppress some of the 
mineralogical fine points, considering type A’s as small, porous, unmelted condensates, and type B’s as large, 
less refractory, and melted since condensation. 
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Figure 5: Cartoon of the hot inner nebula region we adopt for this paper. Region A, between 
r, and r A , is hot enough to evaporate all silicates and refractory oxides. Region B, between 
ta and r B , is hot enough to evaporate all ferromagnesian silicates and metals, but will allow 
refractory oxides to remain solid. Nebula models indicate r A ~ 0.2 AU and r# ~ 2 AU in 
the TO. 5 year timeframe, and decrease quickly afterwards. The nebula evolves inwards at t 
CAIs drift relative to the gas at Vd, and meter-sized particles, with surface mass density 
drift at the larger speed V L . The gas has surface mass density a g and is turbulent, with 
diffusion coefficient V. Silicate “dust” represents f s u of the total disk mass density. 

(CAIs) to evaporate. Beyond that radius, to some other radius r B , CAI minerals can exist 
as solids, but normal ferromagnesian silicates (chondrules and matrix) cannot. In this region 
the temperature is nearly constant - buffered by condensation of ferromagnesian silicates at 
the coolest vertical levels (Bell et al. 1997, Cassen 2001). As the nebula cools, the radial 
boundaries of these regions sweep inwards. While many earlier models seemed to indicate 
that the inner nebula was never hot enough to evaporate silicates (see discussion in Wood 
and Morfill 1989), more recent models seem to allow this. For instance, model calculations 
by Ruden and Pollack (1991), Bell et al. (1997), Stepinski (1998), and Cassen (1996, 2001) 
indicate that ta ~ 0.3 AU and r B ~ 2 — 3 AU at 10 4 years, decreasing to r A ~ 0.01 AU 
and r B ~ 0.3 AU by 10 5 years or so. Note especially that Bell et al. (1997) show that 
temperatures are higher for lower a values, at any given mass influx rate. The very fact that 
the vast bulk of meteoritic materials are isotopically normal to a very high degree suggests 
initial vaporization of a significant fraction of their precursors (eg., Kerridge 1993). We will 
adopt these constraints in our model of the CAI source region. Formation of CAIs only very 
near to the sun (eg. Shu et al. 1996) would restrict the source to something like 0.1 AU or 
less. In our model calculations, we will treat two idealized particle size ranges: “Type A”, or 
< 1mm diameter objects and “Type B”, or > 5mm diameter objects - while larger examples 
of both CAI types are known to exist in smaller numbers. We tentatively associate “Type 
A” condensates with source region A and “Type B” igneous CAIs with source region B. 

2.2.1 Enhancement of silicates and carbon in the source regions 

A new aspect of the model we advance here is a large enhancement in silicate abundance 
in the inner solar system resulting from evaporation of rapidly drifting meter-sized agglom- 
erates of very primitive solid material (which accreted further out, at lower temperatures). 
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This short-circuit for transporting silicate material into the inner solar system determines the 
local mass fraction of vapor in region A from which CAIs eventually form, and/or the abun- 
dance of unvaporized CAIs in region B (both relative to nebula hydrogen). This drift-based 
enhancement process has a limited duration - once the inner nebula cools below the evapo- 
ration temperature of silicates, drifting meter-sized particles are lost directly into the sun. 
Thus, this enhancement has a natural association with the early formation ages associated 
with CAIs. 

The subject of the accumulation of macroscopic particles and planetesimals is a very 
complicated one, and a full description of this process is beyond the scope of this paper. 
The reader is referred to Weidenschilling (1989, 1997, 2000, 2002), and Weidenschilling and 
Cuzzi (1993), for descriptions and some typical recent results. While specific results are 
highly dependent on model assumptions, some generalizations can be attempted. 

For the purpose of this paper, we draw the following generalizations from previous studies: 

(1) The formation of meter-sized particles is relatively easy (in ~ 10 3 — 10 4 years at several 
AU) in either turbulent or non- turbulent nebulae, but growth beyond meter size can be 
severely hindered in turbulent environments. For instance, equation 12# of Weidenschilling 
(1989) leads to a maximum particle radius of ~100 cm for a — 10 -4 , by substituting V t = 
V g = ccc 1 / 2 (the disruption strength is, of course, quite uncertain). In fact, since planetesimals 
form rapidly in current (nonturbulent) models (eg. Weidenshilling 2000), some level of 
turbulence might be required merely to delay the process by the 10 6 year times observed. 

(2) For weak but finite nebula turblence such as we assume, the radial drift velocity of 
meter-sized particles is not decreased by collective effects because the midplane density is 
never large. (3) The particle size distribution for meter-and-smaller size particles can be 
approximated by a powerlaw of the form dN(r ) — r -3 d(logr); this implies roughly equal 
mass per logarithmic radius interval up to meter-sized particles - that is, over 5-7 decades of 
particle size. 2 Based on the width and form of the rubble particle size distribution mentioned 
above, we estimate that during the early CAI-forming stage, meter-sized particles represent 
approximately 10% of the local mass density wherever silicates can be solid. 

Drifting meter-sized silicate agglomerates, rather than being “lost into the sun” during 
this stage, evaporate as they traverse these inner hot regions. To estimate how much mass 
is lost by evaporation, we calculated the evaporation rate of pure forsterite using vapor 
pressure curves of Mysen and Kushiro (1988) and Nagahara et al. (1994), and expressions 
from Supulver and Lin (2000) and Lunine et al. (1991). The evaporation time of a meter- 
radius boulder in our nominal nebula (at pressure of about 1 mbar) varies from 4600 years 
at 1550K to 50 years at 1700K (see also figure 10 of Nagahara and Ozawa 1996). The time 
for such a boulder to drift into the sun from rg = 1 AU is about 100 years. The drifting 
particles have a collisional lifetime that is less than' this, and their fragments will drift more 
slowly while evaporating more quickly. Also, forsterite is the most refractory ferromagnesian 
silicate, so realistic material will evaporate faster. For the purpose of this paper we assume 
the drifting boulders all evaporate and their constituent vapors are homogenized in regions 
A and B. More detailed calculations of how the silicate vapor source function varies with 

2 The particular powerlaw seen in Weidenschilling’s various model results is probably the result of a 
model assumption regarding ejecta size distributions, and the relative fraction of meter-sized particles could 
be higher or lower. 
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radius might be appropriate for incorporation in a more detailed numerical model. 

Because the inward loss rate of vaporized silicate material due to nebula viscous evolution 
into the sun is so much smaller then the arrival rate in large particles, mass conservation 
causes a large enhancement relative to solar abundance of silicate-forming material in either 
the vapor phase, or as small residual refractory objects. This is the radial analog of a 
scheme advocated by Cassen (2001) in which vertical sedimentation played a similar role in 
enhancing heavy element abundance. A similar effect was discussed by Supulver and Lin 
(2000) and Cyr et al. (1998) for drifting water ice particles at a larger radial distance. Below 
we present a toy model which assesses the significance of this process for the CAI factory. 

Consider the hot inner solar system as sketched in figure 5. The region between rg 
and ta is hot enough (>1300K) to evaporate ferromagnesian silicates (95% of rock-forming 
materials; Larimer 1989), but not refractive silicates or oxides - which survive as particulates 
or, perhaps, molten droplets under certain conditions. Silicates enter the region at r B . 
Primordial material (eg., heterogeneous, cometary-like, perhaps even interstellar-like IDPs) 
making its way into this region for the first time probably contains no CAI minerals as such, 
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B. However, because region B is simultaneously cool enough for refractories to condense, 
a haze of fine CAI-type mineral grains might perhaps appear quickly, and without regard 
for the standard condensation sequence (J. Paque, personal communication, 2002). This 
population might be related to “fine-grained” CAIs, and of course would also be subject to 
subsequent melting processes of as-yet unknown nature. The innermost nebula between r* 
and r A is hot enough (1900K) to evaporate everything which reaches it. Vaporized material 
drains into the sun at fj, along with the nebula hydrogen. At values of V characterizing our 
models, we can assume that the concentration C of refractory material (relative to the gas) 

- whether small solid particles or vapor - is quickly homogenized throughout the relatively 
narrow regions .4 and B. 

We obtain the CAI material concentration C, as a mass fraction relative to nebula hy- 
drogen, by solving a simple “leaky barrel” mass balance equation for regions A and B as a 
unit, with both large and small particles entering at rg, small particles and vapor diffusing 
outwards from the region at rg, and vapor draining into the sun at r*. 


where 


d /n ; — dm + 

—(2irra g (r, t)ArC) = — 


dm 

dt 


dm + 
dt 


— 27T r g(C 0 a g B(V n B + Vrfg) + /catO'gVgg), 


( 10 ) 

( 11 ) 


riCa gi V ni + r B Va gB — j. (12) 

In these equations, a g , erg, V n , and V d are surface gas and large particle mass densities, nebula 
inward evolution velocity, and large particle drift velocity respectively. The overline indicates 
a spatial average over the source region of width A r. C 0 = f ca ifsii is the nominal, solar 
abundance of refractory-forming materials by mass relative to hydrogen. This equation may 
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be rewritten as 


ra g (r,t)Ar dC_ _ / fca^jVLB \ c ( r x o gx V r * + / Pa gB (r B /X) \ 

rB&gBiVnB + VdB)Co dt \ C 0 O gB {V n B + Kb) ) \ C 0 r B <7gB(VnB + Kb) ) 

We have simplified the diffusion term of the above equation by taking dC /dr ~ C/A where A 
is the (variable) radial gradient scale due to diffusion itself, and assuming C » C 0 (which 
will be verified). 

To solve this equation we first define a dimensionless timescale tl = t/t c , where 


ro g {r,t) A r 

tbOqB (Kb + Kb) ' 


(14) 


In the above expression for t c , the ratio ( ra g (r,t)/r B crgB ) is a factor of order unity which 
depends on the nebula model, and we neglect it here as a first treatment. Then from figure 
. 3. .wfi. estimate. t s 1.1. — . 5\.x. .1 fri. wars jL, .jTj.a-v-J?e thought of as. the. time in which normal 
evolutionary drift replenishes, or removes, material in a region of width A r. 

Secondly, in the radial range of interest, the inward mass flux 2irre g V n is essentially 
independent of radius (figure 2), leading to cancellations in the last term of equation 13 (to 
order unity). Then, with the further definition of C' = C/C 0 , equation 13 simplifies to 


dC 

^T = ( l + + 

(15) 

which has the simple solution 


C' = E- (E- l)e _(1+ ^ )t '. 

(16) 

In the equations above, 

(17) 

f ca.i® lY LB 

1 ~ C 0 o gB { Kb + Kb) 

(18) 

^ ... Va gB {r B fX) _ V _r B 
&gB(VnB + K B ) r B A Kb A 

(19) 


The enhancement E increases as time proceeds and A increases, but approaches a limiting 
value on the timescale f c /( 1 + T^). The term T\ is the ratio of mass inflow due to large 
particle drift, to that due to gas advection of grains. The term T 2 is the ratio of outward 
diffusive mass flux to inward advective flux. First consider the term T 2 ~ r B / A. Inspection 
of the Green’s functions in figure 6 (next section) indicates that A very quickly approaches 
r B , so on the timescale t c , ~ 1 and E « (1 + E\)/2. 

If the drifting large particles were to contribute negligibly to the inward mass flux, T\ < < 
1 and C ~ C 0 « f ca ifsii for CAIs, where f s u « .005 is the silicate fraction relative to the 
gas and / cai « 0.05 is the canonical fraction of silicates which are refractory (Larimer 1989). 
However, estimates indicate that the drifting large particles in fact dominate the inward mass 
flux of solids. The ratio / ca i<J L /C 0 (J 9 s = o L j f s ua gB is the mass fraction of solids residing 
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in meter sized particles, which we argued earlier in this section is plausibly on the order of 
10 -1 for particle size distributions such as derived by Weidenschilling (1997, 2000, 2002), 
if growth beyond meter-size is hindered by energetic collisions in turbulence. A generally 
expressed estimate of E can be derived as follows, assuming V db « V nB at r B (figures 3 and 
4) and using equation 3 for V n B : 


_ fcajVL Y±b 
2 4 C 0 OgB V n B 


Uol VlbTb 1 jVktb ,2. ( r J.Y 
4 C 0 a gB V AacH Aa\H ) 



(20) 


For 7 = 2 x 10 -3 and a ~ 3 x 10~ 4 — 3 x 10 -3 , as explored in this paper, J57 ~ 10 — 100. 
This simple estimate is surely uncertain at the factor of several level. It is interesting that 
a silicate-to-hydrogen enhancement factor in this range has been inferred for CAI formation 
from detailed chemical and mineralogical modeling by Alexander (2003). 


2.2.2 The role of carbon: chemical and isotopic properties of the CAI factory 

We cigcchr'i'xOtfc, ' quantati very, ' hovv' thfc fact' that the drifting uiaiorxar is millions of years' 
younger, less processed, and far more primitive than “chondritic” can help explain some of the 
chemical and isotopic properties of CAIs and other early, high-temperature condensates. The 
material is in all likelihood enriched in carbon by an order of magnitude over Cl (Jessberger 
et al. 1988, Jessberger and Kissel 1991, Lawler and Brownlee 1992), and a significant amount 
of it is in phases which remain solid to high temperatures (Wooden 2002, Fomenkova 1997). 
Instead of increasing the nebula oxygen fugacity as silicates evaporate, much of the silicate- 
related oxygen could combine with the carbon to produce CO and C0 2 , in relative abundance 
perhaps 1000:1 (Lodders and Fegley 2002; see their figs. 6 and 10, for enhanced Fe/H). This 
may help explain two different aspects of CAI formation at the same time. 

The first effect is chemical: a carbon sink for oxygen (Connolly et al. 1994) may help 
explain how CAIs can form in a gas of nearly solar oxygen fugacity (Beckett et al. 1988) 
even after evaporation of all this silicate material. The second effect is isotopic. There 
are now several indications that the tendency of CAIs and other early condensates to be 
enriched in 16 0 relative to all other solid objects, in a mass-independent fashion, reflects their 
condensation from an ls O-enriched gas, rather than from 16 0-rich solid precursors (Nittler et 
al. 1997, Hiyagon and Hashimoto 1999, Krot et al. 2002). Several ways of increasing the 16 0 
content of the gas from which CAIs condense, in the presence of symmetric molecules such as 
CO2, were reviewed by Thiemens (1996, 1999). The asymmetry provided by incorporation of 
the heavy isotopes 17 0 and 18 0 into excited state C0 2 or CO3 allows them to radiate away 
their energy faster and remain bound more frequently, thus preferentially removing 17 0 and 
18 O from reactive states in the gas and leaving the gas enriched in 16 O (see also Gao and 
Marcus 2001) 3 . 

We note that this entire set of environmental conditions - enrichment in silicate vapor, 
but with solar oxygen fugacity and available C0 2 to affect the isotopic balance - is due to 
evaporation of drifting material. The lifetime of this environment is limited by the cooling of 

3 Alternate theories exist which operate by elevating the heavy oxygen abundance of the rest of solar 
system solids relative to CAIs (Clayton 2002, Lyons and Young 2003, Young and Lyons 2003); however, this 
perspective must still cope with the apparent absense of appropriately 16 0-rich solid precursors (Nittler et 
al. 1997) 
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the inner nebula below the evaporation temperature of ferromagnesian silicates and refrac- 
tory carbon, which happens in the first 10 5 years of the life of the nebula. These conditions 
are thus only associated with the formation age and environment of CAIs and other old, 
high-temparature condensates. 

Summarizing section 2.2, the CAI source regions could be some combination of an annulus 
at the condensation temperature of these minerals, which is some fraction of the width of 
the inner evaporated zone (region A), and a wider, cooler region (region B) where CAI 
minerals are all that remain of evaporated, meter-sized, silicate particles which drift into 
the region from their accretion regions farther out. The boundaries of both zones move 
inwards with time as the nebula cools over perhaps 10 5 years. These source regions, in 
which CAIs are either produced or “refined”, have elevated fractional abundances of CAI 
material: C/C 0 » 1. In a general way, we will explore the hypothesis that “type A” CAIs 
are formed (by condensation) in source region A, and that “type B” CAIs are melted and 
modified in source region B. 


2.3 Radial diffusion in the nebula 

We use a semi-analytical approach which is simpler, and more amenable to assessing wide 
ranges of parameter space, than a full numerical solution to the coupled viscous evolution and 
diffusion difference equations. Clearly, more detailed numerical models would be desirable. 
Our approach makes use of Green’s functions G(r, r 0 ; At). The Green’s functions we use are 
the formal solutions to the radial diffusion equation in 2D cylindrical geometry, assuming 
an instantaneous, narrow cylindrical source of width A r << r 0 at some radius r a . They 
give the resultant diffusion profile at all radii r after some timestep At. We use a At which 
is short compared to the nebula evolution time, but much longer than would be required 
for numerical iteration solutions (Table 1). It is important to note that the quantity which 
diffuses is the local concentration C and not the local mass density Ca g {r,t) of CAIs and 
their vapor presursors; the equation we solve, and to which the Green’s functions apply, is 
(Stevenson 1990, Bockelee-Moravan et al. 2002). 


dC_ 

dt 


a g r dr 



+ So, 


(21) 


where dC/dt is the substantial derivative, and S 0 is the source term. We use a timestep 
which is short relative to the nebula evolution time, because we also need to iterate radial 
drift of the particles in each timestep. The scheme is as follows: (1) obtain nebula properties 
(a g (r, t ) and V n (r, t)) using the LBP solutions; (2) add the CAI source function S 0 , if active, 
to the local CAI concentration C(r, f); (3) convolve with the appropriate Green’s function 
to calculate the diffused CAI distribution at time t + At] (4) calculate the particle drifts 
relative to the gas over time At, and shift particle mass densities (and thus concentrations) 
accordingly; (5) use the gas radial drift velocities to shift CAI concentrations radially; and 
(6) repeat at next timestep. In spirit, this is similar to the operator-splitting approach of 
Bockelee-Moravan et al. (2002). We assume that the source S 0 is only active for 10 5 years, 
approximately the time it takes the inner nebula to cool below the condensation temperature 
of silicates. We selected our code timestep of At — 5 x 10 4 years at the upper end of the range 
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of t c (the viscous replenishment time of an annulus of width A r; equation 14). Selecting the 
smaller end of the range of t c would produce larger CAI concentrations. 

The form of the Green’s function will depend on the radial dependence of the diffusivity 
V. Carslaw and Jaeger (1948, pages 108-109) derived a Green’s function under the condition 
of radially constant diffusivity 4 : 


-(r 2 +r*)/(4T>At) 


( 'TTp \ 


(?i(r,r 0 ;At) 2 tt r 0 J ° \2VAt) 


( 22 ) 


where I 0 is the modified Bessel function of zero order. The argument of I 0 is normally quite 
large in this application, and numerical difficulties arise in the Bessel function itself. These 
are avoided by making use of the large argument form h{x) « e x /\/2 ttx, which is adequate 
over the parameter range of interest, to combine offsetting terms and rewrite the Green’s 
function as: 


Gi(r, r 0 - At) 


/ rjr 

V AirVAt 


e -(r-r 0 )*/{4VA t) 


(23) 


Equation 23 above is similar to the more familiar ID cartesian version, except for a factor 
accounting for the cylindrical geometry. In our numerical models we use the exact small 
argument form of I a when appropriate. 

For the constant ac 2 case, in which viscosity u T varies as r 3 / 2 , we needed to derive a 
different form of the diffusion solution, one which assumes the proper radial variation of the 
diffusivity. This turns out to be 


G 2 (r,r 0 ; At) = 2i rr 0 



=- 4r D /2 ( rl/2+r o /2 )/ T, r>t / g r l/4 r l/4 3/2' 


irV D t 


- 


'Dpt 


(24) 


where Vp is the diffusion coefficient or viscosity at rp, and h(x) is the modified Bessel 
function of third order. The large argument form of G 2 is: 


G 2 (r, r a ;At) 


(r 0 /r) 7/8 (r D /r 0 ) 3 / 4 f 4r 2 (l - (rjr) 1 / 4 ) 2 ^ 
y/^CDAt 6XP V V(r/r D )^At y 


(25) 


The solution for the diffused concentration at subsequent time t + At is the convolution 
of G(r, r 0 ; At) with the radial distribution of concentration at time t: 

rT-max 

C(r, t + At) = / ( C(r 0 , t ) + S(r 0 , t))G(r, r 0 ; A t)dr 0 (26) 

** r min 


where C(r 0 ,t) is the concentration distribution and S(r 0 ,t) is the source at time t, if any. 

Typical Green’s functions for viscosity laws used to generate figures 1 and 2, for several 
different timesteps, are shown in figure 6. It can be seen that on timescales of 10 4 years, 
the radial gradient scale approaches r for both viscosity laws, as argued in our discussion of 
enhancement of the source region (section 2.2). Application of these Green’s functions over 
a finite radial grid (r min , r max ) leads to some mass loss at each timestep, which we correct for 

4 We have rewritten the source “strength” Q of Carslaw and Jaeger in terms of a source concentration 
C 0 , as Q = 2-nr 0 C 0 Ar. 
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Figure 6: Sample Green’s functions for the two viscosity laws we model (figures 1 and 2), 
with source at r 0 =l AU, for three different timesteps. Solid curves: constant viscosity u T = 
10 15 cm 2 sec _1 (G i); dashed curves: constant ac 2 , or v? = 4 x 10 14 (r/lAU) 3 ^ 2 cm 2 sec -1 (G 2 ). 
The three timesteps are 10 3 , 10 4 , and 10 5 years, with the shortest timestep corresponding to 
the narrowest Green’s functions. 
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by enforcing mass conservation across the diffusion substep only. Mass is lost only by radial 
inward drift of the nebula, and the particles in it, into the sun. As a check on the Green’s 
function approach, Taku Takeuchi (personal communication, 2003) kindly ran his model, 
which solves the diffusion equation using difference equations in the normal way, using a 
cylindrical delta function for a source term; the results were in good agreement with our 
Green’s function solutions when the same viscosity laws were used. 

For this initial exploration of the process, we treat the two potential sources (regions A 
and B) separately. The primary difference in the two source types is their radial extent. 
Region A is between 0.01 and 0.3 AU; region B is between 1 and 2 AU. Both sources are 
assumed to be “active” for a duration of 10 5 years. 

3 Results 

Model calculations of the CAI fractional abundance relative to total silicate abundance are 
presented in figures 7 - 14. In all cases, the source regions are assumed to be active for 
10 5 years. The quantity plotted is the concentration - the fraction of “CAIs” relative to 
the local abundance of silicates. Volume fractions of CAIs are what is actually observed, 
but mass fractions are comparable since CAI and chondrule densities are similar. The local 
abundance of silicates is taken to be a constant fraction f s u =0.005 of the local nebula gas 
density as determined by our LBP solutions (figure 1, equations 6 and 8). This simplification 
neglects the variations in silicate/gas abundance due to the diffusion, drift, and evaporation 
processes we discuss above (which only operate for a relatively short time and long before 
primitive body accumulation actually occurs) and decoupling of ~10% of solids in the m-size 
range. However, future models should account for differential evolution of particles and gas 
on longer timescales (eg., Stepinski and Valageas 1996, 1997, Weidenschilling 2002). •• 

The figures have a common format: mass fractions obtained from the constant viscosity 
models are the solid curves; those from the constant ac 2 models, in which viscosity increases 
as r 3 / 2 , are the dashed curves. Results are plotted every 0.5 Myr, from 0.5 Myr to 3 Myr. 
The final curve in each set, at 3 Myr, is shown in a heavier line weight. At the top of each 
plot is a horizontal bar indicating the radial extent of the main asteroid belt. Results are 
shown for source regions A and B separately, for two different particle sizes (diameters of 1 
and 5 millimeters, but one case is shown for 10 mm diameter), and two different values of 
a: 3 x 10 -4 and 3 x 10 -3 . CAIs with larger sizes than these can be found for both major 
types; however, some of the larger size CAIs may well be of a lower density because of a 
high initial porosity, or have a significant irregularity, so would behave aerodynamically like 
correspondingly smaller particles. Thus we select values close to the mode particle sizes. 
Different enhancement factors E are adopted in the different cases to keep the results in the 
ballpark of the observations: to wit, large, 5mm “Type B” CAIs are found (in CV meteorites) 
at the 6% level by volume; small, 1mm “Type A” CAIs are found in essentially all meteorite 
classes at about 1% by volume (Brearley and Jones 1988, McSween 1977a, b). 5 

Small CAIs: We first consider the smaller sized particles, and for these, we first consider 

5 Brearley and Jones (1988) tabulate all inclusions (their page 3-190) but McSween distinguishes between 
true CAIs and other inclusion types, which are less refractory and might have had a longer lasting, or more 
widespread, ’’factory” of their own; thus we match only the true CAI abundances. 
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Figure 7: CAI model concentrations as a function of nebula radius and time, for small 
CAIs (diameter = 1mm), and for our two viscosity/surface density models (solid: constant 
viscosity; dashed: constant ac 2 . This figure is for a Region A source ( r 0 < 0.3 AU), active 
for 10 5 years. Values are plotted every 0.5 Myr up to 3 Myr (which is plotted with a heavier 
weight). . V v. .. 

the Region A source. Figures 7 and 8 show results for 1mm diameter particles, for a =; 
3 x 10 -4 and 3 xlO -3 , respectively. For the a = 3x 10 -4 case, figure 7 shows that Region A 
is able to supply the observed abundance of small CAIs across the asteroid belt - in fact, well 
into the outer solar system - given an enhancement factor in the range of 100-200, depending 
on viscosity model. This is not an unreasonable value of E for this a (equation 20). This 
wide dispersal of small particles is compatible with the results of Bockelee-Moravan et al. 
(2002). For the constant viscosity model, the volume fraction in small CAIs does not vary 
strongly with time. Thus, meteorites formed across the asteroid belt and over the entire time 
period between 1.5-3 Myr might be expected to have the same type and relative abundance 
of small CAIs. For the constant ac 2 model, the CAI volume fraction decreases more rapidly 
with time, but a smaller enhancement factor might apply depending on the time between 
CAI formation and parent body accretion. Here and in all other cases to be discussed, 
models of intermediate radial dependence of viscosity and surface mass density, naturally, 
would produce intermediate results. 

Figure 8, also for a Region A source and mm-diameter particles, is for a larger a = 
3 x 10 -3 . The constant viscosity model is marginal here even for the rather generous en- 
hancement factor of 500, which itself is difficult to justify for this value of a given the simple 
enhancement model we developed (equation 20). The constant ac 2 model does not appear 
to be a contender for this a. 

Consider now the Region B source, and the same small, mm-diameter particles. Figures 
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Figure 8: CAI model concentrations as a function of nebula radius and time, for small 
CAIs (diameter = 1mm), and for our two viscosity /surface density models (solid: constant 
viscosity; dashed: constant ac. This figure is for a Region A source (r 0 < 0.3 AU). Values are 
plotted every 0.5 Myr up to 3 Myr (which is plotted with a heavier weight). An enhancement 
of 500 is probably unrealistic for this value of a = 3 x 10~ 3 . 
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Figure 9: CAI model concentrations as in figures 7 and 8 for small CAIs (diameter = 1mm), 
This figure is for a Region B source (1 < r 0 < 2 AU). Values are plotted every 0.5 Myr up 
to 3 Myr (which is plotted with a heavier weight). 
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Figure 10: CAI model concentrations as in figure 9, for small CAIs (diameter = 1mm), This 
figure is for a Region B source (1 < r Q < 2 AU). Values are plotted every 0.5 Myr up to 3 
Myr (which is plotted with a heavier weight). 

9 and 10 show the results for a = 3 x 10~ 4 and 3 x 10 -3 respectively. The constant 
-viscosity model can satisfy the observations for either value of a , and with considerably 
smaller values of E. This is primarily because source region B is larger in area than Region 
A. As in the Region A case, the temporal variation of volume fraction is fairly small for the 
constant viscosity case, and is larger for the constant ac 2 case. Given the lower enhancements 
required by a Region B source, even the constant ac 2 case might be able to hit the mark in 
the appropriate time window (with a larger enhancement factor). 

Concerning the smaller, “type A” CAIs overall, while region B can satisfy the observa- 
tions with lower enhancement factors, the plausible physical association of region A with a 
condensation origin of many small CAIs leads us to favor it as the source region for “type 
A” CAIs. Appropriate enhancement factors are somewhat larger than our crude estimates 
(equation 20) but not implausibly so. However, enhancement factors required for the larger 
a value (a = 3 x 10~ 3 ) are not consistent with this a (equation 20), so the lower a is pre- 
ferred (eg. figure 7). Given that small CAIs are fairly widely seen in similar proportions in 
all meteorite classes, a situation closer to the constant viscosity model seems to be preferred 
for either Region A or B. 

Large CAIs: Now we address the larger CAIs, taking them to be of 5 mm diameter 
typically as representative of the Type B CAIs which are only seen in one meteorite class 
(the CV chondrites). The data indicate that these large CAIs constitute around 6% by 
volume in the CV chondrites (and essentially zero in other classes). Again, we first look 
at Region A (0.01-0.3 AU). Figure 11 shows that a Region A source with a = 3 x 10~ 4 
can only match the observed abundances if the enhancement factor (2000) is well above the 
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Figure 11: CAI model concentrations as in figures 7 and 8, for large CAIs (diameter = 5 
mm), This figure is for a Region A source (r 0 < 0.3 AU)., and the required enhancement 
(2000) is probably not plausible. 
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Figure 12: CAI model concentrations for large CAIs (diameter = 5 mm). This figure is for a 
Region A source ( r 0 < 0.3 AU) with a = 3x 10“ 3 ; an enhancement this large would be very 
hard to justify for this a (equation 20), and even at that, the abundances are significantly 
lower than observed in the >lMyr timeframe. 
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Figure 13: CAI model concentrations as in figures 9 and 10, for large CAIs (diameter = 5 
mm). Curves are plotted every 0.5 Myr from 0.5-3 Myr. This figure is for a Region B source 
(1 < r a < 2 AU). The’ observed abundances (6%) may be satisfied by the constant viscosity 
model with an enhancement factor which is quite reasonable for this a. There may even be 
enough latitude in E for the steeper viscosity law to satisfy the observations, as well. 
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Figure 14: CAI model concentrations as in figure 13, for large CAIs (diameter = 5 mm), 
This figure is for a Region B source (1 < r 0 < 2 AU). For this a, enhancements as large as 
1000 are questonable, so neither viscosity law provides a good candidate. 

maximum plausible value (equation 20). Even for this large enhancement, the steeper nebula 
viscosity models (dashed lines) are unable to match the observed volume fractions after 1.8 
Myr or so. The situation is even worse for the larger a of 3 x 10~ 3 (figure 12), where this 
large enhancement factor is even less reasonable. Neither viscosity model can retain 5mm 
particles for 1 Myr in the observed abundances. So, region A is an unlikely source for “type 
B” CAIs within the context of this model. 

On the other hand, the Region B source (1-2 AU) is able to supply the larger particles in 
.the observed abundances because of its larger area and its closer proximity to the asteroid 
belt. Figure 13, for a = 3x 10~ 4 , shows that at lelast the constant viscosity model matches 
the observations with a plausible enhancement. However, even the Region B source is 
questionable if a is as large as 3 x 10 -3 (figure 14), because for this a an enhancement factor 
of 1000 is problematic (equation 20) and even then, after 1.5 Myr the model abundances are 
smaller than the observations of large CAIs require, even for the constant viscosity model. 

Some Type B CAIs are even larger - in the 1 cm diameter range. Unfortunately, no 
quantitative data on CAI size distributions exists in the literature to our knowledge. Figure 
15 shows results for Type B CAIs as large as 1 cm in diameter. The larger particles are 
clearly not so easily retained over periods as long as 3 Myr; however, for a value of E which 
is not implausible for a = 3 x 10~ 4 , some of these objects can be retained until 1.5-2 Myr 
or so. at least by the constant viscosity model. The abundance decreases even more rapidly 
than for 5mm particles, over only a few 10 s years. 

For 5-10mm diameter CAIs, the time variation of the modeled CAI abundances is much 
stronger than the radial variation, and because it decreases montonically with time, a pre- 
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Figure 15: CAI model concentrations as in previous figures, for especially large CAIs (di- 
ameter = 10 mm), This figure is for a Region B source (1 < r„ < 2 AU). There are no 
quantitative estimates for the fractional abundance of cm-diameter objects. 

diction of the model would clearly be that the CV chondrites in which large, “tyP e B” CAIs 
are found should be accumulated somewhat earlier than other chondrite groups, rather than, 
say, at the same time but at a different distance from the sun. 

4 Discussion and Conclusions 

We briefly summarize our results, and mention some additional implications, predictions, 
and tests of the model. 

We have shown that, in a weakly turbulent protoplanetary nebula, outward diffusive mass 
flux across a concentration gradient can offset inward mass flux of l-5mm diameter particles, 
by nebula advection and gas-drag-driven radial drift, for times on the order of 10 6 years. 
The relative abundance of small (1mm diameter) or highly “fluffy” particles (which behave 
aerodynamically like smaller, solid particles) is so slowly changing in the 1-3 Myr timeframe 
that it is not surprising to find them ubiquitously distributed in all meteorite classes. The 
radial distribution of small objects is fairly uniform, and some redistribution even into the 
outer solar system is indicated, consistent with earlier results for even smaller particles by 
Bockelee-Moravan et al. (2002). The model thus offers some interesting possibilities for 
whether or not small CAIs would be observable in comets or IDPs; generally speaking, 
they would be at lower abundances (sometimes far lower) than in meteorites, but for some 
parameter choices, they might be present in observable quantities. This result should be 
treated cautiously, because our model calculates the particle radial drift rate using a constant 
pressure gradient term 7 (see section 1.2), which will not be the case towards the outer edge 
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of the nebula. 

During the very same initial 10 5 years when CAIs seem to have formed, the inner nebula 
was hot enough to evaporate quite a large mass in inwardly drifting meter-size boulders, 
thereby greatly enhancing the abundance of silicate material relative to solar. We present a 
quantitative estimate of the enhancements which might result from this process, and show 
that plausible enhancements provide quantitative agreement with observed CAI abundances 
for a range of turbulent viscosities and assumptions as to the radial extent of the source re- 
gion, as well as with independent estimates based on CAI formation and mineralogy (Alexan- 
der 2003). We suggest that carbon enrichment of these early boulders, relative to Cl, might 
help to explain aspects of meteorite chemistry which argue for solar oxygen fugacity, and for 
mass-independent fractionation of oxygen isotopes. This silicate-and-carbon enhancement 
process ceases when the inner nebula cools off - well before chondrules form and primitive 
bodies accumulate. Thus, the “CAI factory” environment is naturally associated with the 
early, high temperature condensates which show these properties. 

The association of small “type A” condensates with source region A, and large, “type B” , 

’ gn oous CAIs with P oop-rpc to b° plctusiblo. Tb.sit i.s ? r. rojpon A .(ooT>dR T ?SRtio n . «•. 

zone) source is able to explain the observed abundance of 1mm or “fluffy” condensate parti- 
cles with enhancements in a plausible range. Based on the enhancement factors involved, it 
seems that only a region B source is likely to explain the observed abundances of these larger 
CAIs. A slightly cooler region B is also consistent with the somewhat less refractory miner- 
alogy of type B CAIs. Of course, condensate particles from region A would need to diffuse 
outwards through region B, allowing those which remain there for a longer time to have the 
opportunity to coalesce, melt, and/or equilibrate with a somewhat cooler gas. Small “type ■ 
A” condensate CAIs we observe in meteorites today would be those which traversed region 
B quickly, incurring the least processing subsequent to their formation. The larger E values 
associated with region A formation of small CAIs might be suggesting that evaporation of 
drifting rubble culminates primarily at or near the region A boundary, rather than uniformly 
throughout regions A and B. 

5-10mm diameter particles are retained less easily than 1mm diameter particles, and 
their abundance is rapidly decreasing in the 1-3 Myr timeframe. Thus, a prediction and test 
of the scenario and model presented here is that the CV class of meteorites, which alone 
contain the large, 5-10mm CAIs, should have accumulated before the other classes. It is a 
bit of a puzzle as to why it is said that large CAIs are not only rare, but nonexistent in other 
chondrite classes. Even a drop by an order of magnitude or more should leave some rare 
occurrences of large CAIs in later-forming classes. This might indicate a layer of complexity 
not in the current model. 

Overall, the viscosities/diffusivities toward the higher end of our range (a = 3 x 10~ 3 ) 
make the scenario less credible, because the high E required becomes inconsistent with those 
very a values. This is because we assume that viscosity evolves the nebula as well as diffusing 
the particles. Larger a values decrease the nebula gas density more quickly, increasing 
particle stopping times and inward radial drift speeds faster than outward diffusion can 
proceed. A model in which the diffusivity and nebula evolution were decoupled might not 
incur this difficulty, and indeed, turbulent diffusivity is always larger than turbulent viscosity 
(Prinn 1990). 

The meteorite record is still incompletely sampled, and some specifics of current gener- 
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alizations may be vulnerable. For instance, Kimura et al. (2002) have observed abundant 
CAIs in a group of related ordinary (H) chondrites - in spite of the current generalization that 
CAIs are rare in ordinary chondrites. These appear to be the small, type A CAIs which are 
widely distributed in carbonaceous chondrite classes - but in larger relative abundance than 
normally found in ordinary chondrites. Moreover, the matrix material in this H chondrite is 
like that of C03 chondrites - introducing the possibility that the nature of the fine-grained 
material, which is expected to diffuse the most intimately with the gas, might be a more 
important factor in CAI abundance than the chondrule type with which it is mixed. This 
leaves a diffusive mixing scenario still attractive, but hints at spatial and/or temporal vari- 
ations beyond the scope of this paper. Clearly, more observations of diverse meteorite types 
•will be enlightening. 

Our model is certainly amenable to improvement in many ways. Intermediate, but still 
analytical, forms of nebula viscosity and mass density distributions can readily be explored 
with our Green’s function approach, but significant improvement beyond that will quickly 
bring a number of complex parameter variations into play. For instance, a realistic numerical 
diffusion equation, with coupled accretion and drift, would be required to explore radially 
and temporally variable viscosity and diffusivity. At the same time, a radially variable 
pressure gradient berm (7) could be incorporated to address the transport of material to the 
outer solar system. Another key process which should be explored is the incorporation of a 
“loss” term into the model (eg., Cassen 1996, 2001), whereby grains can be incorporated into 
much larger particles as they evolve outwards, thereby ceasing their outward evolution and 
probably settling to the midplane as part of the large particle. Perhaps already solid, mrri- 
to-cm-size CAIs are less efficiently accumulated than vaporized ferromagnesian material, 
which quickly freezes out into tiny grains or onto the surfaces of already existing grains 
and particles as the enhanced silicate vapor inside region B diffuses outwards across its 
condensation boundary. The possibilities for recycling of silicate material through the hot 
inner zone are intriguing. Finally, realistic timescales for mineralogical and chemical changes 
by reaction with nebula gas could be combined with this model to assess the timescales on 
which particles evolve radially through regions of different temperature. This would result 
in another set of testable predictions. Clearly, there is much work remaining to be done. 


28 



Acknowledgements: We are grateful to Robbins Bell, Pat Cassen, Harold Connolly, 
Tom Greene, Larry Grossman, Sasha Krot, Glenn MacPherson, Bjorn Mysen, Katharina 
Lodders, Mark Marley, Kevin McKeegan, Julie Paque, Mark Thiemens, John Wasson, Diane 
Wooden, and Kevin Zahnle for helpful discussions concerning CAIs, carbon, isotopic frac- 
tionation, silicate evaporation, and nebula properties. We thank Robert Last and Robert 
Hogan for computational and graphics support, Pat Cassen for careful reading of an earlier 
version and several useful suggestions, and Ignacio Mosqueira and Steve Desch for other 
helpful comments. We thank Taku Takeuchi for running a numerical test of our Green’s 
function approach. This research was supported by Grants to JNC from the Planetary Ge- 
ology and Geophysics Program and the Origins of Solar Systems Program. This research 
has made ample use of the Astrophysics Data System Query and Abstract service. 



Table I 

Model Parameters 


Parameter Definition 


M 0 

M d 

K 

c 

H 

a 

V, iat 

tmax 

A t 
tb 

r A 

n 

r D 


Stellar mass 

Protoplanetary nebula mass (t = 0) 

Solar system specific angular momentum 

Gas molecule thermal speed 

Nebula vertical scale height 

Turbulent difFusivity parameter 

nebula turbulent difFusivity and viscosity 

Duration of hot inner nebula CAI source region 

Code evolutionary timestep « t c 

inner boundary for normal silicates (t < t max ) 


„,-.r .. - 


(x ^ 4 \ 

V'- ''max j 


inner boundary of nebula 
angular momentum or centrifugal radius 


Value or range 
2x“i0 33 "g 
0.2 M 0 

9.5 x 10 19 cm 2 sec“ 1 
1 km sec"" 1 

0.05 x radius 

3 x 10“ 4 - 3 x 10~ 3 

acH 

10 5 yr 

5 x 10 4 yr 

2.0 AU 

no ATT 

v> . ^ Ao 

0.01 AU 

4.5 AU 


Other quantities 


Parameter 

Definition 

Reference 

?*p, dp 

particle radius and diameter 



t. 

particle stopping time due to gas drag 

eqn. 

1 

. fail , /cot 

silicate/hydrogen, CAI/silicate mass fractions 

0.005, 0.05 

C, C 0 

CAI material concentration in general, and nominal 

eqn. 

14 

E ''' ' • 

Enhancement of CAI concentration relative to nominal 

eqn. 

17 

St 

particle Stokes number 

eqn. 

Iff 

v d , V dB . 

CAI particle radial drift velocity in general and at r B 

. eqn. 

2,3 

Vn, V nB 

nebula radial evolution velocity in general and at r B 

eqn. 

2,3 

v K ,n 

Keplerian velocity and orbit frequency 



v L ,v LB 

large particle drift velocity in general and at r B 



7 

pressure gradient term affecting gas drag 

eqn. 

2 

(7, j ^2, &ga O gB 

nebula gas surface densities 

eqns 6, 8 

OL,Olb 

large particle surface density in general and at r B 



FuK 

parameters determining enhancement E of C 

eqn. 

18,19 
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